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We have developed a new theoretical formalism for phonon transport in nanostructures using the 
nonequilibrium phonon Green's function technique and have applied it to thermal conduction in 
defective carbon nanotubes. The universal quantization of low-temperature thermal conductance in 
carbon nanotubes can be observed even in the presence of local structural defects such as vacancies 
and Stone- Wales defects, since the long wavelength acoustic phonons are not scattered by local de- 
fects. At room temperature, however, thermal conductance is critically affected by defect scattering 
since incident phonons are scattered by localized phonons around the defects. We find a remarkable 
change from quantum to classical features for the thermal transport through defective CNTs with 
increasing temperature. 

PACS numbers: 44.10.-|-i, 61.46.Fg, 63.20.Mt 



The discovery of the quantization of phonon thermal 
conductance as well as quantization of electrical con- 
ductance has had a great impact on mesoscopic and 
nanoscopic physics 0,0, Q • Since the first observation 
of quantized thermal conductance kq = ir'^k^T/Sh in 
dielectric mesoscopic wires by Schwab et al. signif- 
icant effort has been expended on exploring new ma- 
terials exhibiting kq. Carbon nanotubes (CNTs) are 
expected to be potential candidates to measure quan- 
tized thermal conductance |^ because they arc ideal 
one-dimensional phonon conductors with a large phonon- 
mean-free path H, S 111- Recent sophisticated exper- 
iments have in fact succeeded in measuring quantized 
thermal conductance in CNTs 0. 

Recently, concern has been raised that various intrin- 
sic features of pure CNTs arc lost because of the pres- 
ence of defects in synthesized CNTs. According to recent 
molecular-dynamics (MD) simulations, the thermal con- 
ductivity in CNTs at room temperature decreases dra- 
matically with increasing defect density UlI ll^ . 
Similar to room-temperature thermal transport in defec- 
tive CNTs, low-temperature thermal transport in CNTs 
is predicted to be markedly affected by defects, and quan- 
tized thermal conductance might be destroyed due to the 
presence of the defects. To the best of our knowledge, the 
influence of defect scattering on low-temperature thermal 
transport through defective CNTs has not been studied 
thus far. The central aim of this Letter is to cl arify the 
effects of structural defects, such as vacancies hj and 
Stone- Wales (SW) defects [l5|, on quantized thermal con- 
ductance in CNTs. Contact scattering resistance is not 
considered here to focus on the intrinsic conductance of 
defective CNTs. 

In order to clarify the effects of structural defects, a 
reliable theory of thermal transport through a nanostruc- 
ture is needed. Rego and Kirczenow heuristically derived 



a novel formula for thermal current through a mesoscopic 
dielectric wire placed between hot and cold heat baths [jj] : 
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where fL{R) (w) is the Bosc-Einstcin distribution func- 
tion of equilibrium phonons with energy fko in the left 
(right) lead with temperature (Tfl), and Ci'-^) is the 
phonon transmission function. The Landauer formalism 
has some remaining issues to be resolved. First, it is 
not applicable to an interacting phonon system includ- 
ing phonon-phonon interactions. A second issue is that 
there is no systematic procedure to obtain Ci'-^) for ac- 
tual nanomatcrials with complex structures on an atomic 
scale. Thus far, (^{u) has been determined using elastic 
models for classical sound waves. Of course, such a treat- 
ment is not suitable for nanoscale systems with complex 
atomic structures. In this Letter, we develop a new for- 
malism overcoming these remaining issues by utilizing 
the nonequilibrium Green's function (NEGF) technique. 
The advantage of the NEGF formalism is that local phys- 
ical quantities, such as the nonequilibrium phonon den- 
sity and the local thermal current, can be calculated. 
This Letter is a first report of the application of the 
NEGF formalism to the study of the influence of defect 
scattering on phonon transport in CNTs. 

The system Hamiltonian J^/fgys is described as the sum 
of the harmonic term 
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and the anharmonic term J^anh.- Here, Sia{t) is an op- 
erator in the Heisenberg picture for atomic displacement 
from equilibrium along the a direction of the ith atom 
with mass Mi. pia{t) is a momentum operator conju- 
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gated to the displacement operator Sia(t), and kia,j[3 rep- 
resents the spring constant between the ith atom in the 
a direction and the jth atom in the f3 direction. The sys- 
tem Hamiltonian is assumed to be divided into five parts: 
J^L+J^LS + '^s + '^RS + '^R- Here, Jf^/u. is the 



Hamiltonian for the left/right thermal lead, =^ is that 
for the scattering region, and J^ls(RS) is the Hamilto- 
nian for the coupling between the scattering region and 
the left (right) lead. The anharmonic term J^^nh is also 
assumed to exist only in the scattering term J^s ■ Differ- 
ent temperatures, Tl and Tr(< Tl), are assigned to the 
left and right regions of the system, respectively. 

The thermal current flowing through the interface be- 
tween the left lead and the scattering region can be cal- 
culated from the time evolution of the energy of the left 
lead: Jth = — (-^) = ■^([■^, ^ys]), where the bracket 
(• • • ) denotes the nonequilibrium statistical average of 
physical observable. The thermal current Jth is rewrit- 
ten as 
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using the greater and lesser Green's functions associ- 
ated with the contact between the left lead and the 
scattering region: i?iG> ^^(t, t') = {sia{t)sji3{t')) and 
ihGf^ -p{t,t') = {sjf3{t')sia{t)). Since the Green's func- 
tions depend only on the time difference in a steady 
state, it is convenient to work in Fourier space (w space). 
Therefore, in the steady state, the thermal current is ex- 
pressed as 
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where the relation Gf^ ^^{—uj) 



Gff3,io,{'^) is used. 



In accordance with the similar procedure for the NEGF 
formalism for electronic transport [3, Eq. O can be 
straightforwardly rewritten as 
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where the boldface quantities represent matrices with ba- 
sis in the scattering region. This is a general expression 
for the thermal current beyond the Landauer formula 



to the left /right lead and is calculated using the mode 
matching method |l7l | modified for phonon transport. 

Phonon-phonon scattering is not important for CNTs 
below 300 K, which is of interest here 0, 0, @l- the 
case of coherent phonon transport, the lesser and greater 
Green's functions satisfy the Keldysh equation: 



(7) 



The retarded/advanced Green's function G^J°'{lj) for the 
scattering region satisfies the Dyson equation: 



■/a 



(8) 



where D is the dynamical matrix derived from the second 
derivative of the total energy with respect to the atom 
coordinates in the scattering region, and M is a diag- 
onal matrix with elements corresponding to the masses 
of the constituent atoms. In this work, the total energy 
of the CNTs is determined from the Brenner bond-order 
potential [l8| . 

Substituting Eqs. © and (0) into Eq. lO, the thermal 
current in Eq. (|3J) is reduced to the Landauer formula 
in Eq. The phonon transmission function C,{lo) in 

Eq. 1^ is also expressed explicitly as 



C(c.) = Tr [Tl{lo)G's{oj)Tr{lo)G%{uj)] 



(9) 



This expression of C,{lo) is equivalent to that derived by 
Mingo and Yang in a different way In the limit 

of the small temperature difference between the hot and 
cold heat reservoirs, Tl - Tr <^ T = {Tl + Tr)/2, the 
thermal conductance k = Jth/ (Ti — Tr) is given by 
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If the phonon transmission is perfect for all acoustic 
modes in the low-temperature limit T ^ 0, Eq. H10() 
is given as a form of an elementary integration that can 
be integrated analytically, and the thermal conductance 
is quantized as Mkq = M {n'^ k'^T / 3h) , where M is the 
number of acoustic modes. 

We now apply the formalism developed here to the 
thermal transport in defective CNTs. Figuren]represents 
the phonon transmission function for the (8,8)-CNT with 
and without defects. The dashed curve is the transmis- 



for coherent phonon transport. In Eq. ©, G|(tj) is the sion function Cp(^) for the perfect (8,8)-CNT without 
greater/lesser Green's function for the scattering region, any defects, and it displays a clear stepwise structure that 
and S<^j^(a;) is the greater/lesser self-energy due to cou- gives the number of phonon channels. In the low-energy 
pling to the left /right lead, which is given by region below 2.4 meV, being the energy gap of the low- 

est optical modes, the dashed curve shows Cp(^'^) = 4 
indicating the number of acoustic branches correspond- 
ing to longitudinal, twisting, and doubly degenerated 
fiexural modes. Reflecting the perfect transmission for 
all acoustic modes, the thermal conductance k{T) shows 
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(uj) = ^iUl/rH + 1/2 ± 1/2)Tl/r{uj). (6) 



where T,''/^/r{lo) 
is the retarded/advanced self-energy due to the coupling 



Here, Tl/r{uj) = i[S2/^(w) 
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FIG. 1: (color online). Phonon transmission function (^{hui) 
for the (8,8)-CNTs. Red and blue curves represent 
for the vacancy and SW defect, respectively. The dashed 
curves are for the perfect CNT. The top-left inset describes 
the phonon density around the SW defect at 7.0 meV indi- 
cated by the blue arrow, and the bottom-right inset is that 
around the vacancy at 11.6 meV indicated by the red arrow. 
The red shading on the atom spheres indicates the phonon 
density. The red shading is not on the same scale for the two 
insets. (The figures of insets are dropped from this preprint 
due to large sizes.) 

four universal quanta. Akq — A{n'^ k'^T / 3h) , in the low- 
temperature limit (see Fig.lSJ. 

The transmission function Cvac{fii^) for the (8,8)-CNT 
with the vacancy and Csw(fi^) for the (8,8)-CNT with 
the SW defect are described by the red and blue curves 
in Fig. n] respectively. The Cvac/sw(^'-^) is dramatically 
deformed from (p{huj) owing to defect scattering partic- 
ularly at high energies. However, it remains unchanged 
in the low-energy region. This is because the long wave- 
length acoustic phonons with low energies in the CNTs 
are not scattered by local defects such as vacancies and 
SW defects. This leads to the important conclusion that 
the thermal conductance of CNTs exhibits Akq at cryo- 
static temperatures even when the CNTs include a small 
amount of defects (see Fig.|21l. 

The transmission function Cvac(fi^) shows some dips at 
particular energies, which are clearly distinguished from 
the dips originating from the van Hove singularity of the 
optical phonon branches. These dip positions coincide 
with the peaks in the local density of states (LDOS) 
around the vacancy (not shown). The appearance of 
LDOS peaks means that the phonon density is highly lo- 
calized around the vacancy. The transmission dips arise 
from the scattering of incident phonons from the lead by 
the phonon localized states. The bottom-right inset in 
Fig. n] shows the phonon density around the vacancy at 
11.6 meV indicated by the red arrow in Fig.^ The LDOS 
peak at 11.6 meV is located at the lowest position and 



FIG. 2: (color online). Low-temperature thermal conduc- 
tances in (6,6)-, (8,8)- and (10,10)-CNTs. Red and blue 
curves represent k{T) / hq for the CNTs with the vacancy and 
SW defect, respectively. The dashed curves are for the perfect 
CNT. 

largest intensity among the peaks. Similarly, Csw(^) 
has some dips at particular energies coinciding with po- 
sitions of the LDOS peak due to the SW defect. The dip 
at 7.0 meV indicated by the blue arrow lies at the low- 
est position among the LDOS peaks associated with the 
SW defect. As shown in the top- left inset in Fig. ^ the 
phonon density at 7.0 meV is strongly localized around 
the SW defect'. 

Substituting the C(^'^) obtained into Eq. l(TTi|l . we can 
determine the thermal conductance k{T) as a function of 
temperature T. The low-temperature behavior of k{T) 
normalized to the universal quantum kq for (6,6)-, (8,8)- 
and (10,10)-CNTs with and without defects are shown in 
Fig El The red (blue) curves represent Kvac(sw)/'^o foi' 
the CNTs with the vacancy (SW defect). The dashed 
curves represent Kp/^o for perfect CNTs. As discussed 
above, all thermal conductance curves approach 4 in the 
limit of T ^ 0, even in the presence of local defects. 
In the temperature region shown in Fig. |21 the thermal 
conductance ksw/kq is slightly lower than Kvac/Ko- This 
is because propagating phonons in this region are domi- 
nantly scattered by the phonon localized state associated 
with the SW defects at 7.0 meV, and not by the phonon 
localized states associated with a vacancy at the higher 
energy of 11.6 meV (see Fig^. 

We next describe the CNT-diameter dependence of 
defect scattering on thermal conductance for moderate 
temperatures up to 300 K. Figure O shows the ratios 
Kvac/Kp and Ksw/np for (6,6)-, (8,8)- and (10,10)-CNTs 
as a function of T. The black-, red-, and blue-solid 
(dashed) curves are the ratios Kvac(sw)/'*p for (6,6)-, 
(8,8)- and (10,10)-CNTs with the vacancy (SW defect), 
respectively. All the Kvac(SW) / '^p curves decrease rapidly 
with increasing temperature and become nearly indepcn- 
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FIG. 3: (color online). The temperature dependence of the 
ratio Kvac(sw) /'i^p for the CNT with the vacancy (SW de- 
fect). The black-, red-, and blue-solid (dashed) curves are 
Kvac(sw)/'^p for (6,6)-, (8,8)- and (10,10)-CNTs, respectively. 



dent of the temperature at '^300 K. The terminal value 
of 'vvac(sw)/'^p s-lso decreases as the CNT gets thinner. 
In other words, the influence of defect scattering in thin 
CNTs on the thermal conductance is more significant 
than that in thick CNTs. Interestingly, Kvac/^p is clearly 
lower than Ksw/'tp at moderate temperatures, which is 
in sharp contrast with the low temperature case. That is, 
the incident phonons are scattered more strongly by the 
vacancy than by the SW defect. This result is consistent 
with previous results of classical MD simulations at room 
temperature |l3l |. Therefore, the most important finding 
in this study is that defect scattering, which is respon- 
sible for thermal transport, changes remarkably from a 
quantum to a classical nature with increasing tempera- 
ture. 

Finally, we note the electron contribution to the ther- 
mal conductance in metallic CNTs. The ballistic elec- 
trons in metallic CNTs without any defects contribute 
four quanta (= 4ko) to the thermal conductance |^|^, 
according to the Wiedemann- Franz relation between elec- 
trical conductance and electronic thermal conductance. 
Even in the presence of the SW defect in metallic CNTs, 
electronic thermal conductance is quantized as 4ko since 
the SW defect does not modify the electrical conductance 
near the Fermi level [20I lilj. On the other hand, with 
the vacancy, the electronic thermal conductance is re- 
duced because the conduction electrons arc scattered by 
the dangling cr-bond states around the vacancy with an 
energy close to the Fermi level j2r)| . 

The infiuence of defect scattering on phonon transport 
in CNTs differs from that on electron transport. For 
example, as discussed above, the propagating electrons 
in CNTs with a vacancy are strongly scattered by the 
dangling cr-bond states |20j , whereas the long wavelength 
phonons arc not perturbed by the vacancy. 



In conclusion, we have developed a new formalism for 
phonon thermal transport in nanostructures using the 
NEGF technique. Applying this to CNTs with local 
structural defects, a vacancy and a SW defect, we found 
for the first time that a remarkable change in defect 
scattering, i.e., defect-dependent thermal conductance of 
CNTs from a quantum to a classical feature, occurs with 
increasing temperature. Our formalism opens the way 
for a complete understanding of the underlying physics 
of phonon transport in nanostructures under various con- 
ditions. Finally, we discussed the advantages of our for- 
malism based on the NEGF over other theories. Our the- 
ory, in a straightforward manner based on the Feynmann 
diagram technique |2^ . takes into account the phonon- 
phonon scattering effect that is found to play an impor- 
tant role in thermal tr ansp ort above room temperature 
in a recent experiment 
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